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INTRODUCTION 

Supersonic  combustion  in  turbulent  hydrogen-air  systems  is  a  topic  of  current  interest 
because  of  its  relevance  to  applications  such  as  the  development  of  one-stage-to-orbit 
vehicles  that  employ  air-breathing  propulsion.  Previous  research  in  the  present  program 
established  that  in  applications  of  this  type  the  combustion  often  occurs  in  the  reaction-sheet 
regime.  Current  work  has  therefore  been  focused  on  determining  the  structures  and 
combustion  limits  of  reaction  sheets  in  these  hydrogen-air  systems.  Influences  of  different 
descriptions  of  the  chemical  kinetics  of  combustion  are  being  investigated  with  the  objective 
of  identifying  the  simplest  kinetic  scheme  that  gives  results  that  are  sufficiently  accurate  for 
use  in  engineering  computations  for  the  combustors  of  applied  interest.  In  addition,  the 
compressibility  effects  are  being  studied  in  various  geometrical  configurations. 
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RESEARCH  OBJECTIVES 


The  objective  of  this  research  is  to  improve  understanding  of  the  chemical  kinetics 
and  fluid  dynamics  of  turbulent  combustion  in  high-speed  flows.  Supersonic  combustion 
in  hydrogen-air  mixtures  is  being  addressed  by  theoretical  approaches  that  distinguish 
between  reaction-sheet  and  distributed-reaction  regimes.  The  work  seeks  to  identify  effects 
of  compressibility  in  turbulent  combustion,  methods  for  including  compressibility  in 
theoretical  analyses,  and  reduced  chemical-kinetic  mechanisms  appropriate  for  supersonic 
combustion.  The  results  may  help  to  enhance  capabilities  of  reasonable  computations  of 
high-speed  turbulent  reacting  flows. 

ACCOMPLISHMENTS 


Boron  Combustion 

Although  boron  combustion  is  not  part  of  the  present  program,  the  AFOSR  grant 
preceding  the  present  one  concerned  boron  combustion,  and  publication  of  results  from  that 
work  is  being  completed  under  the  present  grant  Items  1-3  of  the  list  of  publications 
concern  articles  on  results  of  the  research  on  boron  combustion. 


In  the  most  recent  work  on  this  subject  attention  was  focused  on  including  low- 
temperature  oxidation  of  boron  with  a  liquid  oxide  present  on  the  particle.  Although 
diffusion  of  B  through  the  liquid  layer  was  established  to  be  most  important  at  higher 
temperatures,  detailed  consideration  of  Russian  experimental  results  strongly  suggested  that 
at  lower  temperatures  the  diffusion  of  O2  becomes  more  important.  An  expression  was 
derived  for  the  temperature  at  which  this  change-over  in  the  controlling  mechanism  occurs. 
The  principal  accomplishments  of  our  previous  work  can  be  summarized  as  follows: 

A  description  of  ignition  and  combustion  of  boron  particles  was  developed  with  the 
aim  of  obtaining  a  theory  that  is  as  simple  as  possible  yet  consistent  with  available  data. 
For  the  ignition  stage  the  model  involves  equilibrium  reactive  dissolution  of  B  in  the  thin 
8203(1)  layer  to  form  dissolved  BO,  surface  attack  of  BO  by  02(g)  and  by  H20(g)  to  form 
802(g),  and  HOBO(g),  respectively,  vaporization  of  B2O3,  and  reaction  of  H2O  with 
8203(1)  to  form  HOBO(g).  For  the  combustion  stage  the  model  involves  clean-surface 
attack  of  B  by  02(g)  and  B203(g)  followed  by  a  gas-phase  diffusion  flame  between  B  and 
O2.  Values  of  relevant  rate  constants  were  identified  that  achieved  agreement  with 
experiments.  Early  laser-ignition  experiments,  not  previously  explained  quantitatively, 
were  employed  for  obtaining  needed  dry-gas  rate  parameters.  New  experiments  on  ignition 
and  combustion  of  boron  suspensions  in  hot  combustion  products  of  a  flat-flame  burner 
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were  performed  to  test  adopted  values  of  wet-gas  rate  parameters.  Results  provide  a  basis 
for  calculation  of  ignition  and  combustion  of  boron  in  propulsion  applications. 

Turbulent  Combustion  Regimes  in  Supersonic  Combustion 

In  previous  work,  calculations  had  been  completed  of  Damkolhler  and  Reynolds 
numbers  for  the  anticipated  range  of  supersonic  combustion  conditions  showing  turbulence 
Reynolds  numbers  to  range  from  103  to  107  and  Damkdhler  numbers  from  10  to  104.  The 
latter  results  motivated  the  present  focus  on  the  reaction-sheet  regime. 

Counterflow  Hvdrogen-Air  Flame  Structure 

The  structure  and  extinction  of  counterflow  diffusion  flames  of  hydrogen  and  air 
were  investigated  for  pressures  from  0.5  to  10  atmospheres  and  for  initial  temperatures 
from  300  K  to  1200  K  as  described  in  items  4  and  5  of  the  publications.  Numerical 
integrations  were  performed  for  air-side  strain  rates  from  60  s-1  to  extinction.  The 
numerical  results  were  compared  with  predictions  of  an  asymptotic  analysis  that  involved 
reduction  to  one-step  chemistry  through  introduction  of  steady-state  and  partial-equilibrium 
approximations.  Reasonable  agreement  was  found  for  concentrations  in  the  main  reaction 
zone  at  low  strain  rates. 

When  these  flame-structure  computations  were  first  completed  and  extinction 
predictions  were  compared  with  available  experimental  results  (from  the  NASA  Langley 
facility),  it  was  found  that  the  predicted  strain  rate  at  extinction  exceeded  the  measured 
value  by  nearly  a  factor  of  two.  However,  our  computational  results  motivated  refined 
experiments  at  Langley  that  produced  the  revised  experimental  extinction  strain  rates,  for 
different  hydrogen  mole  fractions  Xh2  in  hydrogen-nitrogen  fuel  mixtures,  as  listed  in 
Table  1.  Also  shown  in  Table  1  are  predictions  of  two  other  computational  groups.  It  is 
seen  in  Table  1  that  our  predictions  are  in  excellent  agreement  with  experiment  and  that  our 
agreement  is  better  than  that  achieved  by  the  other  groups.  We  attribute  this  success  to  our 
careful  selection  of  the  best  elementary  rate  constants  available  in  the  current  literature.  It 
seems  especially  noteworthy  that  for  the  hydrogen-air  system  our  prediction  preceded 
experiment  and  was  subsequently  confirmed. 
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Table  1:  Comparison  between  experimental  and  one-dimensional  numerical 
evaluations  of  extinction  strain  rates  (s-1)  at  atmospheric  pressure 
with  the  feed  streams  at  300  K. 
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In  the  application  of  flamelet  models  to  turbulent  flames  it  may  be  valid  and  certainly 
efficient  to  consider  simplified  descriptions  of  both  the  chemical  kinetic  and  the  fluid 
mechanical  processes.  Therefore  as  an  effort  complementary  to  that  described  earlier  an 
analysis  of  hydrogen-air  flames  in  counterflowing  configurations  with  reduced  chemical 
kinetic  and  simple  transport  descriptions  is  underway.  Programs  permitting  the  limiting 
cases  of  frozen  and  equilibrium  as  characterized  by  the  flame  sheet  approximation  are 
completed  and  incorporation  of  a  reduced  mechanism  to  deal  with  finite  rate  chemical 
effects  is  presently  underway.  The  results  of  these  simplified  descriptions  will  be 
compared  with  the  full  calculations  discussed  earlier  in  order  to  assess  the  possible  errors  in 
applications  to  turbulent  flames.  In  this  same  regard  attention  is  being  devoted  to  the 
special  features  associated  with  the  incorporation  of  flamelet  models  into  turbulent  flames  in 
high  speed  flows  in  particular  with  the  influence  of  fluctuations  of  temperature  and  pressure 
in  the  two  streams.  A  review  which  includes  some  comments  on  these  studies  is  to  appear 
soon  as  publication  5. 

In  a  more  recent  study  we  have  employed  asymptotic  methods  (publication  6)  to 
analyze  the  displacement  effects  produced  by  heat  release  in  counterflow  flames.  This 
work  addresses  the  significant  corrections  called  for  to  interpret  properly  experimental 
results  in  such  flames,  in  particular  on  extinction  strain  rates.  Introduction  of  these 
corrections  was  shown  to  bring  theory  and  experiment  into  good  agreement 

In  even  more  recent  work  we  have  included  nitrogen  chemistry  in  our  computations 
because  we  learned  at  the  previous  AFOSR  contractor's  meeting  that  concerns  existed 
about  the  possibility  of  nitrogen  chemistry  significantly  influencing  flame  structure  and 
extinction.  Our  recent  results  demonstrate  these  influences  to  be  negligible.  For  example, 
the  effect  on  the  maximum  temperature,  which  we  find  to  be  greatest  when  we  add  10~3 
ppm  of  NO  to  a  hot  air  stream,  is  always  less  than  a  30  K  change  at  a  temperature  in  the 
range  of  2800  K.  For  these  reasons  we  do  not  currently  plan  to  study  effects  of  nitrogen 
chemistry  further.  Instead,  our  current  direction  seeks  greater  simplification  in  kinetics. 

Extinction  computations  were  performed  with  all  H2O2  chemistry  removed  from  the 
system,  and  the  influence  was  found  to  be  negligible.  Recently  computation  schemes  with 
four-step  and  with  three-step  reduced  mechanisms  have  been  implemented  in  our 
calculations.  Preliminary  results  indicated  that  both  of  these  reduced  mechanisms  yield 
excellent  agreement  with  results  of  computations  based  on  the  full  mechanism.  Work  is 
proceeding  to  test  three-step  and  two-step  approximations  for  achieving  maximal 
simplicity.  We  already  know  that  a  one-step  scheme  will  not  give  good  extinction  results. 

Compressibility  Effects 

In  previous  work  on  this  project  a  theoretical  analysis  of  the  in  viscid  flow  between  a 
porous  plate  and  a  parallel  impermeable  plate  was  performed  for  small  values  of  the  ratio  of 
the  plate  separation  distance  to  the  lateral  extent  of  the  plates,  for  both  planar  and 
axi symmetric  geometries,  as  described  in  item  7  of  the  publications.  The  problem  of 
computing  the  flow  field  was  reduced  to  the  solution  of  a  single  integral  equation,  which 
was  accomplished  numerically.  The  ratio  of  specific  heats  y  was  a  parameter  of  the 
solution,  and  parametric  results  were  obtained  from  y  =  1.0  to  y  =  1.67.  The  flow 
exhibited  choking  at  a  critical  value  of  the  lateral  extent  of  the  plate,  in  the  vicinity  of  which 
the  Mach  number  approaches  unity.  The  results  are  needed  in  providing  external 
boundary-layer  conditions  for  studying  the  flame  structure  in  the  viscous  region  between 
two  counterflowing  streams  when  compressibility  is  important. 
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It  was  observed  that  the  same  formulation  that  was  employed  in  this  study  also  can  be 
applied  to  a  model  for  the  flow  field  in  a  laterally  burning  rocket  motor.  For  this  reason, 
such  an  extension  was  completed,  although  it  is  somewhat  peripheral  to  our  central 
objective.  The  results  are  presented  in  publication  7.  Future  work  on  compressibility 
effects  is  intended  to  return  to  counterflow-like  configurations,  considering  structures  that 
may  arise  in  vortices. 

SUMMARY 

These  continuing  studies  are  helping  to  improve  our  understanding  of  structures  of 
reaction  zones  in  high-speed  flows.  Effects  on  descriptions  of  turbulent  combustion 
remain  under  active  consideration.  In  addition,  investigations  of  compressibility  effects 
and  of  reduced  chemical-kinetic  mechanisms  continue  to  be  pursued  to  improve 
understanding  of  turbulent  supersonic  combustion. 
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